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ABSTRACT 

The boundary layer for a swirling flow of incompressible 
fluid in a convergent cone has been investigated. Theoretical 
calculation with boundary conditions more appropriate to 
physically existent situations discount the existence of 
"super velocities" within the boundary layer. Parallel 
experimental investigations demonstrate an interdependence of 
core and boundary layer flows which precludes the maintenance 
of initial "ideal flow" conditions. 
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NOMENCLATURE 



velocity in x direction - axial velocity 
velocity in 6 direction - tangential velocity 
velocity in z direction 
distance along the wall 
distance perpendicular to wall 

initial ratio of tangential to axial velocities 

free stream velocity in x direction 

derivative with to X of U 

free stream tangential velocity 

derivative with respect to X of V 

acceleration vector 

velocity vector 

vorticity 

viscosity 

kinematic viscosity 
density 



shear at the wall 
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I. 



INTRODUCTION 



The determination of velocity profiles in the laminar 
boundary layer of a convergent nozzle has been approached 
analytically in several different manners. Taylor [1] con- 
sidered the case of a predominant tangential flow with a 
secondary axial flow. In his analysis, Taylor assumed that 
the tangential flow exceeded the axial flow to the extent 
that the axial component in the inviscid core could be 
neglected. After applying the Pohlhausen method to the 
mementum-integral equations and obtaining solutions for various 
axial positions, Taylor concluded that a condition could 
arise when most of the outflow from the nozzle was fed by a 
boundary-layer current close to the wall of the swirl chamber. 

The radial pressure gradient, which held the particles in their 
circular path, acted upon the retarded boundary layer to drive 
the fluid along the surface to the exit. 

Wilks [2] developed a similar approach that led to the 
possibility of a 'super-velocity' existing in the boundary 
layer of a convergent conical nozzle. This 'super-velocity' 
was a velocity within the boundary layer that exceeded the 
free-stream velocity. Wilks considered the flow in the cone 
to be that of an inviscid core with a thin boundary layer. 

The inviscid core had two components of flow - a uniform 
axial flow and a free vortex with its center along the center 
line of the nozzle. 

In this paper an analogous approach was utilized. How- 
ever, the effect of the boundary-layer growth upon the tangential 
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and axial velocities in the free stream flow was considered. 
Subsequently, no super-volocities were found to exist. Further- 
more, experimental results prove that full scale effects 
predominate in this problem rather than just boundary-layer 
effects . 
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II. ANALYSES 



A. WILKS THEORY 

Using the Wilks approach, a curvilinear coordinate system 
was defined with the x-axis along the wall and the z-axis 
perpendicular to the wall of the nozzle (Figure in Appendix 
A) . The shape of the nozzle was determined by the surface of 
revolution of an arbitrary curve r(x). Wilks developed an 
approximate method of solution using the momentum equations 
(Appendix A) 
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and the corresponding integrated forms of the momentum 
equations 
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V = V(x) = b—} ^ = 0; b-JL = 0; when z ->■ <» 

r (x) 3z 3 z2 

The distances in the above equations were non-dimensionalized 
by defining n = z / 6 and X = x / c , where 6 was the local 
boundary-layer thickness and c was the distance along the 
generator to the nozzle to its apex. 

To obtain a solution of the momentum-integral equations, 
it was necessary to assume a velocity profile within the 
boundary layer. A profile that would also satisfy the energy- 
integral equations should produce some improvement over a 
solution to the momentum-integral equation alone. Wieghardt 
[3] developed such a method with the use of two free parameters 
and an eleventh-degree polynomial. Applying the free para- 
meters and the polynomial to the u-velocity profile and ass\am- 
ing that the v-velocity profile did not change significantly 
with X, the following profiles were obtained: 

u/U = f]^(n) + L 2 *f 2 (ri) + L 2 *f 3 (n) (8) 



where 

f^(n) = 1 - (1 - n)®(l + 8n + 36n^ + 120n) 

f 2 (n) = (1 - n)®n(l + 8n + 36n2) 

f_(n) = -(1 - n)®n2(l + 8n) 

3 
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and 



v/V = g(n) = 2n - 2n^ + n'^ (9) 

An evaluation of the axial velocity profile at the wall leads 
to physical connotations for the shape factors, and L 2 ; i.e. 

includes the wall shear stress and L 2 involves the axial 
pressure gradient. 

T _ Tq 6 _ 66 dP 

1 yU 2 2c yU ' dX^ 



Using the coefficients Al, A2, A3 and A4 (ratios of character- 
istic thicknesses) as defined in APPENDIX A and introducing 
a term 



G = (^) (|) = {^) 



the momentum- integral equations became 
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As a method of simplifying the solution, Wilks assumed that a 
basically linear relationship existed between the momentum 
thickness ('S2x^ boundary-layer thickness (6) and 

between the mixed momentum thickness (^2xy^ boundary- 

layer thickness . The above equations were then combined 
with the definition of G and its derivative. The following 
relationships were obtained 
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A2 = G ^(1 + Al) + A4 - G (g)^ A3 
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In the preceeding analysis, no restriction was made upon 
the shape of the convergent nozzle. The flow conditions in 
a cone, as well as many other shaped nozzles, do not meet the 
requirements for a similarity solution of the related boundary 
layer equations. However, in these cases where the profiles 
vary with axial position in a non-linear manner, this present 
approximate method of solution can be universally used. Wilks 
showed that a solution of the similarity equation could be 
developed for any one point along the surface of a cone, but 
that the coefficients of the Falkner-Skan equation would be 
different for each point. For each point, he noted that a 
close agreement existed between the profiles of the approximate 
method of solution and the profiles of the similarity solution. 

For the particular case of a cone, a linear relationship 
existed between the radius and X. The axial velocity as a 
functions of X was then formulated from continuity, and the 
variations in the tangential velocity were derived from con- 
servation of angular momentum. (The details of these functions 
are given in APPENDIX A.) 
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B. 



NUMERICAL SOLUTION OF WILKS THEORY 



A solution of the momentum- integral equations would lead 
to the values of the shape factors and thus to the axial 
velocity profile. A step-by-step procedure was used to solve 
the governing equations at various values of X. Wilks 
considered the simplified case where the initial boundary- 
layer thickness on the cone was zero. With this assumption, 
G(0) became zero, which led to L2(0) = 0. Equation 13 then 
became a quadratic in L^^CO). The solution for other values 
of X continued in the following iterative manner: 

(1) Using the values of L^CO) and L2(0), the 
thickness coefficients were calculated. 

(2) The velocities were then evaluated. 

(3) Using equation 13, L^^ was calculated. 

(4) Applying equation 12, the value of G* was 
determined . 

(5) Using the first two terms of a Taylor series 
expansion, G(X+dX) = G(X) + G' (X) *dX. 

(6) L 2 was calculated from equation 14. 

(7) The values of the thickness coefficients were 
recalculated with the new values of L^^ and L 2 
and the procedure from step 2 onward was 
repeated for new values of X. 

In a physical system, a boundary layer could easily exist 
before the flow entered the conical region. The momentum- 
integral equations would then need to be solved for a non-zero 
initial boundary-layer thickness. This was accomplished by a 
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slight modification of the Wilks solution. Using the solution 
of the quadratic in L^CO) for zero initial boundary-layer 
thickness as a predicted value of Lj^(O), an iterative method 
would correct this value for the presence of an initial 
boundary-layer thickness. The method pursued was the following 

(1) An initial value for Lj^(O) was assumed. 

(2) The velocities were calculated. 

(3) G and thus L 2 were determined from the value 
of the initial boundary-layer thickness. 

(4) Using the current values for Lj^(O) and L2(0), 
the thickness coefficients were calculated. 

(5) The value of Lj^(O) was determined according 
to equation 13. 

(6) The calculated value of L-j^ was compared with 

the assvied value of If the two values 

were not sufficiently close, a value equal to 
the average of the two was assumed and the 
procedure from step 4 onward was repeated. 

If the result was sufficiently accurate, 
and L 2 were calculated for other values of X 
using the same methods as for a zero initial 
boundary-layer thickness . 

Using the above numerical methods, the momentum-integral 
equations were solved for values of K (swirl factor) of one 
and two. The resultant profiles were compared with the Wilks 
profiles to insure the correctness of the numerical method 
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employed. The solution was then extended for values of K up 
to six with various initial thicknesses and initial velocities 
Several profiles that were obtained in this manner are shown 
in Figures 1-7 . The solutions indicated that the degree of 
super-velocity increased dramatically with K and X. The 
value of the momentum thickness became negative under many 
circumstances, and for high degrees of super-velocity, the 
displacement thickness also became negative. Under all of 
these conditions the boundary-layer growth continued to 
increase with X. The numerical analysis deteriorated for 
values of X much greater than 0.3 in most of these solutions. 
There was also a failure in the numerical solution for ratios 
of tangential velocity to axial velocity in the range of three 
It was believed that these problems were caused by the numer- 
ical method employed and were not inherent problems of the 
analytical model. The quantity (A2 - A4) could be very close 
to zero. Since one of the terms in the calculation of 
required division by the above quantity, overflow in the 
computer could easily be caused when this quantity approached 
zero. A different final arrangement of the momentum-integral 
equations would avoid this problem and allow for a solution 
at higher values of X. 
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FIGURE 1. VARIATIONS IN PROFILE WITH K. X=0 . 1 G(0)=0 

(a) K=0 , (b) K=1 , (c) K=2 , (d) K=3 , (e) K=4 
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FIGURE 2. VARIATIONS IN PROFILE WITH K. X=0 . 1 G(0)=0.52 
(a) K=1 , (b) K=2 , (c) K=3 , (d) K=4 
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FIGURE 3. VARIATIONS IN PROFILE WITH K. X=0.2 G(0)=0.52 
(a) K=1 , (b) K=2 , (c) K=4 
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FIGURE 4. VARIATIONS IN PROFILE WITH X. K=1 G(0)=0.52 
(a) X=0.00 , (b) X=0.05 , (c) X=0.10, (d) X=0.15 

(e) X=0.20 , (f) X=0.25 
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FIGURE 5. VARIATIONS IN PROFILE WITH X. K=2 G(0)=0.52 
(a) X=0.00 , (b) X=0.05 , (c) X=0.10 , (d) X=0.15 

(e) X=0.20 , Cf) X=0.25 
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FIGURE 6. VARIATIONS IN PROFILE WITH X. K=4 G(0)=0.52 

(a) X=.00 , (b) X=.05 , (c) X=.10 , (d) X=.15 , (e) X=.20 
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FIGURE 7. VARIATIONS IN PROFILE WITH G. K=2 X=0.2 
(a) 6{0)=0.00 , (b) 6(0)=0.10 , (c) 6(0)=0.15 
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C. PRESENT THEORY 

The Wilks approach to the momentum equation assumed that 
the tangential velocity at the outer edge of the boundary 
layer could be described as A/r. This was only appropriate 
if the boundary-layer thickness was very small compared to 
the value of r. The n\americal solution of the Wilks method 
led to a boundary-layer growth that could be of the same 
order as r for flow in air where r^ was about one inch and the 
initial axial velocity was one foot per second. In these 
cases, the radius minus the boundary-layer thickness should 
be used to determine the tangential velocity at the outer 
edge of the axial boundary layer rather than using the radius 
alone. Thus V = A/(r - 6) . 

The introduction of this new value of V into the boundary 
conditions of the momentum equations led to a different form 
of the v momentum-integral equation. The derivation of the 
equation is shown in APPENDIX B. This approach required the 
definition of a new thickness coefficient A5 (APPENDIX B) . 

The V momentum-integral equation then became 




In the above equation, Wilks assumption of a basically linear 
relationship between the characteristic thicknesses was used. 
The terms in the braces are added to the Wilks v momentum- 
integral equation by the new assumption, vix., V = A/(r - 6). 

If the boundary layer is thick enough to effect the 
boundary condition for the tangential flow, the calculation 
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of the axial velocity should also take this boundary-layer 
growth into consideration. In using the conditions of conti- 
nuity of flow, the radius minus the displacement thickness 
should be used rather than the radius alone. The displacement 
thickness, by its definition, is the distance the potential 
flow has been displaced by the presence of the boundary layer. 
The effect of this consideration would be to increase the 
axial core velocity as the displacement thickness grew. 

D. NUMERICAL SOLUTION OF PRESENT THEORY 

In order to calculate U and its derivative, U', the value 
of the displacement thickness and its derivative were needed. 

The ratio of the displacement thickness to the boundary-layer 
thickness was one of the integrals used in the calculation 
of the thickness coefficients . The boundary-layer thickness 
at any point could easily be calculated from the value of 
G once the kinematic viscosity and length of the nozzle to 
its apex had been specified. The calculation of the derivative 
of the displacement thickness was not as direct. Wilks based 
his solution on the assumption that a basically linear relation- 
ship existed between the momentum thickness and the boundary- 
layer thickness. An extension of this assumption would be 
to presume a basically linear relationship also existed between 
the displacement thickness and the boundary-layer thickness. 

Thus the ratio of the displacement thickness to its derivative 
should be equal to the ratio of the boundary-layer thickness 
to its derivative. 
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with the new equation for the v momentum integral, the 
numerical solution proceeded in the following manner: 

(1) Using the initial boundary-layer thickness, 
the values of G(0) and L2(0) were calculated. 

(2) A predicted value for Lj^(O) was assumed from 
the Wilks solution for zero initial boundary 
layer thickness. 

(3) The distance along a flat plate that would be 
required to produce a boundary-layer thickness 
of the same size as the initial boundary-layer 
thickness was calculated. Using this distance, 
the derivative of the boundary-layer thickness 
on a flat plate at this point was calculated. 
From this estimate of the derivative and 
assuming an initial value for D1 from the Wilks 
solution, the initial value of the derivative 
of U was evaluated. 

(4) Using the current values of and L 2 / the 
values of the thickness coefficients were 
calculated. 

(5) The value of 6'/8 was calculated, where 



6 ' ,A4 U' _ r' A5 

T~ “ U r “ r/6 - 1' 




(1 + 



A5 



r/6 - 1' 



(16) 
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(6) The new value of was calculated from 



L, = D2-G* r — (2 + Al) + (1 - xi a3) + —1 

^ Lu r u2 6 J 

(17) 

(7) The calculated value of L^^CO) was compared with 
the assumed value. If the two values were not 
sufficiently close, the thickness coefficients 
and velocities were re-evaluated until the 
values of Lj^(O) for tv/o consecutive iterations 
were sufficiently close. 

(8) The value of G' was determined from 



G' = G 




(9) The values of G and L 2 for the next point v;ere 
then calculated. 

(10) The step-by-step evaluation of the momentiam- 
integral equations was continued for new values 
of X by returning to step 4 and remembering 
that step 7 is used only when X=0. 

U and V as functions of X are given in APPENDIX B. 

Using this numerical technique, the momentum-integral 
equations were solved for various values of initial boundary- 
layer thicknesses, initial velocities, and K. Since this 
method required the division by G in the calculation of equa- 
tion 16, the analysis was not valid for the case where the 
initial boundary-layer thickness was zero. However, for very 
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small values of G other than zero, there was no restriction 
on the initial boundary-layer thickness . The solution of 
these equations showed no sign of super-velocity. In several 
cases the step-by-step evaluation continued smoothly for values 
of X up to 0.6. The boundary-layer growth reached a maximum 
value when X was in the range of 0.25 and then decreased with 
X. This numerical analysis appeared to lose its validity when 
the momentum thickness became negative, which occurred at a 
much higher value of X then in the Wilks solution. At this 
value of X where the momentum thickness became negative, dis- 
continuities appeared in the boundary-layer growth. Figures 8-12 
show a comparison between the profiles of the Wilks solution, 
as modified for initial boundary-layer thickness, and the 
profiles obtained by the above method. 

E. DISCUSSION OF THEORETICAL RESULTS 

Wilks based his solution of the flow in the boundary layer 
of a convergent nozzle on several assumptions. He first con- 
sidered the inviscid core as a free vortex with uniform axial 
flox. It was then ass;amed that the boundary layer would be 
thin enough to prevent any effect of its growth upon the core 
velocities . Thus the core velocities were dependent only upon 
the geometric configuration of the convergent nozzle. Wilks 
also assxmed that a linear relationship existed between the 
various characteristic thicknesses of the boundary layer. With 
these assximptions , the solution of the momentum-integral 
equations led to the existence of a super-velocity in the 
boundary layer. 
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For low velocity air flow (required to maintain the 
laminar conditions) the boundary-layer growth v/as not negli- 
gible. By including this growth in the calculations of the 
boundary conditions, the solution no longer led to super- 
velocities. Increased swirl only decreased the ratio of the 
displacement thickness to the boundary-layer thickness 
slightly . 
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G(0)=0.52 6 (0) /r^=0 . 130 (a) Wilks theory , (b) present 

theory 
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FIGURE 9. COMPARISON OF THEORETICAL CURVES. K=1 X=0.2 G=0.52 
6 (0) /r^=0 . 130 (a) Wilks theory , (b) present theory 
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FIGURE 10. COMPARISON OF THEORETICAL CURVES. K=2 X=0.1 
G(0)=0.52 6 (0) /rQ=0 . 130 (a) Wilks theory , (b) present 

theory 
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FIGURE 11. COMPARISON OF THEORETICAL CURVES. K=2 X=0.2 
G(0)=0.52 6 (0)/r^=0.130 (a) Wilks theory , (b) present 

theory 
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FIGURE 12. COMPARISON OF THEORETICAL CURVES. K=4 X=0.1 
G(0)=0.26 6 (0) /rQ=0 . 065 (a) Wilks theory , (b) present 
theory 



35 



III. EXPERIMENTATION 



A. DESIGN OF APPARATUS 

The analytical approach assumed that the flow entering 
the convergent nozzle was a combination of a free line-vortex 
and uniform axial flow. To test the theory, an apparatus was 
required to approach this assumed velocity distribution. There 
are several methods of introducing a fluid into a chamber and 
causing a swirling motion. However, this tangential velocity 
would have to approximate the conditions of a free vortex, that 
is the tangential velocity should be inversely proportional to 
the radius. Several methods were considered, including rotating 
a section of the swirl-producing chamber, drawing the fluid into 
the chamber from a stagnation condition while using vanes to 
create the tangential motion, and rotating blades inside the 
chamber. The approach chosen was to introduce the fluid through 
one tangential port with variable cross-section. In tests of a 
Hilsch tube. Lay [4] used a single circular port to introduce 
a jet into the chamber. This was the only source of fluid in 
his tests, and he reported the velocity profiles that were 
produced in this manner. However, he did not indicate the 
radial profile of the velocities in a region close to the wall. 
The results indicated that the tangential flow might closely 
approximate the conditions for a free vortex over most of the 
radius, departing from these conditions in the center portion. 

The distribution of the axial velocity in Lay's tests suggested 
that the superposition of pipe flow upon this flow could result 
in a nearly uniform axial flow. 
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The source of fluid for this apparatus was a compressed 
air tank. To insure laminar pipe flow of the air at speeds 
of at least two feet }>er second, a pipe diameter of two inches 
was chosen. Axial pipe-flow could be produced prior to entry 
of the axial flow into the swirl chamber. The shape of this 
axial flow could be varied by adjusting the axial position of 
a honeycomb f low-straightener . If profiles different from 
pipe flow were required to achieve the uniform axial profile 
at the entrance of the cone, the flow could be adjusted by 
only passing air through selected annular portions of the 
honeycomb flow-straightener . 

The tangential flow was introduced through a rectangular 
slot in the wall of the swirl geneaator. One wall of this 
slot was tangent to the inner surface of the chamber. The slot 
was 0.25 inches high and the width could be varied from 0.3 
to 1.25 inches. This variable cross-section permitted an 
adjustment in the ratio of the velocity to the quantity of mass 
introduced by the tangential port. This flexibility in the 
system was added to insure that the proper profiles and values 
of K could be obtained. Tests were then performed to determine 
the axial distance required for the development of the appro- 
priate flow conditions. These preliminary tests indicated that 
a range between one and two feet was required for the development 
of the tangential jet into a free vortex. 

Lay's tests indicated that there would be a region close 
to the wall where the tangential flow would be retarded by 
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viscosity. An axial combination of pipe flow and the tangen- 
tially introduced flow would also have a region close to the 
wall where the conditions for uniform axial flow would not be 
met due to the viscous forces present. With these two boundary- 
layer growths under consideration, the outer annular portion 
of the flow was bled off, with only the "flat" inner portion 
of the axial flow entering the cone. This bleed-off was 
accomplished by placing a cone of initial radius 0.77 inches 
at the end of the tube which had an inside radius of one 
inch. The position of the cone relative to the end of the 
tube could be varied in order to alter the amount of flow bled 
off into the atmosphere. The entry region into the cone 
consisted of a one-inch straight section. The purpose of this 
section was to insure the continuation of the proper flow 
characteristics after bleed-off. (See Figure 25) 

The cone, a piece of machined plexiglass, featured a 
semi-angle of 4.4 degrees. The distance along the generator 
of the cone to the apex, termed c in the analysis, was ten 
inches. The truncated cone used in the test apparatus extended 
over 0.35 c. Thus, tests could easily be performed over the 
region where X varied from 0.0 to 0.3. Two hot-wire probes 
were introduced from the exhaust end of the truncated cone. 

One configuration vi?as sensitive to axial and radial components 
alone . The other probe could detect tangential and radial 
velocities. It was assumed that the radial component of flow 
was negligible, especially close to the wall. The traversing 
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mechanism moved in two directions, allowing for motions in 
both the x-direction and the z direction. 

B. EXPERIMENTAL RESULTS 

The development of flow conditions that would approximate 
the desired conditions was basically a trial-and-error pro- 
cedure. By adjusting the slot width, the development length 
for the tangential flow (by changing the length of pipe between 
the swirl generator and the nozzle), the bleed-off level, and 
the radial distribution of the axially introduced flow, it was 
possible to achieve several flow conditions that approximated 
the analytical boundary conditions at the entrance of the cone. 
However, in all the cases tested, the velocity distribution in 
the inviscid core chan ^ed drastically with X. The radial dis- 
tribution of the axial and tangential velocities for several 
values of X with swirl ratios (K) of 0.6, 0.8, and 0.95 are 
shown in Figures 13-24. These profiles indicate several 
different flow relationships. For the case where K=0.6, the 
degree of swirl had very little effect upon the axial flow as 
it moved through the cone. The uniform axial core was gradually 
replaced by a flow that resembled the distribution in pipe flow, 
with most of the flow occuring in the central portion, and the 
perimeter retarded by viscous forces. With the slight increase 
to K=0.8, there was a great change in the radial distribution 
of the axial flow. In this case, the free vortex seemed to be 
maintained for all test stations. At X=0.3, this tangential 
motion had thrown a great deal of the fluid into the peripheral 
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portion of the cone, leading to a high axial velocity near the 
wall. However, this high velocity was still in the core 
( 2 = 0.2 inches) and not in the boundary layer. In the tests 
where K=0.95, the distance over which the tangential flow 
appeared as a free vortex was less then in the previous cases. 
This smaller distance resulted in several changes in the 
tangential velocity with X. As a result of this changing 
tangential profile, there were several in the axial flow. In 
sections where the tangential flow appeared as a forced vortex, 
there was one form of axial flow, where high speed occurred 
away from the wall, while in annular sections where the tan- 
gential flow resembled a free vortex, the high speed axial flow 
occurred nearer the wall. 

C. DISCUSSION OF EXPERIMENTAL RESULTS 

In both of the analytical approaches, the inviscid core 
was considered a combination of uniform axial flow and a free 
vortex. In the test cases, the axial flow at the entrance of 
the nozzle was within five per cent of uniform axial flow. The 
tangential flow approximated a free vortex over most of the 
radius in the entrance region. However, this flow could not 
be maintained down the cone. Swirl induced a high velocity 
axial flow near the wall when the tangential flow was in the 
form of a free vortex. This region of high axial-velocity flow 
was far enough away from the wall (a distance greater than 0.15 
inches) to call it part of the inviscid core. Thus swirl does 
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indeed cause changes in the axial flow, but these changes occur 
in the core and not in the boundary layer as Wilks suggested. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



Previous work in this area considered the effect of 
swirling flow upon the velocity profiles in the boundary 
layer. In this paper, the effect of the boundary layer growth 
upon the core velocities was utilized in the analysis. How- 
ever, the experimental tests showed that this consideration 
was not enough. The swirl also affected the axial core flow, 
and the resultant changes are much more drastic than any varia- 
tion in the boundary layer flow. Rather than applying the 
Navier-Stokes equations to the boundary layer alone, it would 
be much more appropriate to apply them to the whole flow field. 
In addition, the developed theory for the boundary layer assumed 
that a basically linear relationship existed between the 
various characteristic thicknesses. Since these thicknesses 
are functions of the shape factors, which are in turn functions 
of X, it would seem mandatory to extend this analysis without 
these assumptions of linear relationships. However, in order to 
complete such an analysis, more information needs to be known 
about the effect of swirl upon the core velocities. 
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APPENDIX A 



DERIVATION OF EQUATIONS ACCORDING TO WILKS 
The boundary of the convergent nozzle is formed by sur- 
face of revoltuion of a curve r(x), where x is a distance 
measured along the wall. The distance from the wall is z. 
The potential flow has the two components, V and U as shown 
below 




Acceleration in the orthogonal curvilinear coordinates is 
defined as 

® ^ it = If + ^ (l/2g2) - q X t 

In the i . direction this becomes 
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If = X , U2 = 6 
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For axisymmetric flow, all of the partial derivatives with 
respect to the angular displacement are zero. Since velocities 
in the direction perpendicular to the wall (the z direction) 
may be considered small compared to the other velocities, 
the momentum equation in the z direction may be neglected. By 
an order of magnitude analysis, the boundary-layer equations 
for steady laminar flow with constant properties become 
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36 
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I- (r-u) + I- (r-w) = 0 

oX 3 Z 



with the boundary conditions 
u=v=0 on z=0 

u U(x) and v -»■ V(x) = A/r (x) as z 
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with the above boxindary conditions, the pressure gradient in 
the x-direction can then be written 

- i !£ = u — - £l 

p 9x dx r 

Integrating the u-inoraentum equation across the boundary layer 
gives 






9u V dr 
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„ + v2 

r dx “ dx ^ r j 
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from continuity 

w = - i [ I— (r‘u) dz 
r J 0 9x 
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The second term 
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can be integrated by parts . 
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00 ,00 
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The u momentum- integral equation then becomes 

[ fu ^ (r*u) + H ^ (r*u) - U — ^ dz 

J o\ ax r 9x r 9x dx ^ 
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The first integral can be written 
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The integral form of the v-momentum equation is 

O'* ^ s) = - - (If)„ 



r • u) dz 
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Introducing the continuity equation, this becomes 
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Performing an integration of the second term by parts 
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The left hand side of the momentxam-integral equation may now 
be written as 

^ fIf - 

= f“rf a(vu) _ V 1^^ - V iii + V IH - u li 

JoLv^^ ” 9x; 9x 9x r 9x 

+ XIX dz 

r dxj 

0 9 , F Ti 9u V»u 9r , 2’u*v dr \ ^ 

_ S -5F - — ^ -F- aFj 



Since 
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the equation can be written 
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Thus the V momentum-integral equation is 
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with the boundary conditions 
V = 0 
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The u momentum-integral equation is 
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In order to obtain a solution to the momentum-integral 
equations using the approximate method of von Karman, it is 
necessary to assume suitable velocity profiles. A direct 
analogy with the Pohlhausen method would lead to 



where n = z/6 (x) . The continuity equation leads to a stream 

function 

\j; = r«U*h(x,n) 

The axial velocity can then be written 



If the tangential velocity is of the form v = V*g(x/n) , 
the momenttim equations may be written 

h"" + a-h-h" + b*(l-h'2) + d* (l-g2) = O 

g" + a*h*g' = 0 

where a, b, and d are constants as defined by Wilks. Boundary 
conditions require that g" evaluates at the wall is zero. A 
differentiated form of the v-momentum equation gives 
g"' + a* (h' *g' + h*g") = 0 

Since both h' and g" evaluated at the wall are zero, it is 
apparent the g"' evaluated at n = 0 must also be zero. This 
would require the second and third derivatives of g to be 



u/U = f(n,Li) v/V = g(n,L 2 ) 



u = U h(x,n) 
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independent of x. With this understanding, the two parameters 
can be used to obtain the axial velocity profile while the 
tangential profile is considered dependent upon n alone. 

Thus 

u/U = f(n,L^,L 2 ) and v/V = g (n) 

Introducing a dimensionless length X=x/c (c being the length 
of the cone along the generator) , the u momentvim-integral 
equation becomes 
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The V momentum-integral becomes 
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with the boundary conditions 
g = 0 

g 
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* on n = 0 
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Here it was assumed that the boundary conditions could be 
transferred from infinity to the boundary-layer thickness, 
that is where n=l, without appreciable error. 

In an attempt to obtain an improvement in the basic 
approximate method of solution, a profile can be obtained that 
will also satisfy the requirements of the energy-integral 
equation. Weighardt [3] developed such a method using two- 
parameter profiles. This profile can be represented with the 
aid of an eleventh order polynomial. Thus 



u/U = fi(n) + Li*f 2 (n) + L 2 *f 3 (n) 
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Evaluating the functions at the wall (n=0) yields 





thus 
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Also 
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The integrals in the momentum integral equations can be 
expressed as various characteristic thicknesses. 
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With the definition of a new thickness 

,1 



= 6 



lo 



(1-g^) dn 



the momentum-integral equations can now be written 
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Differentiating the appropriate terms, these equations 
simplify to 
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An integration of the polynomials in the above integrals led 
to 

D1 = 0.3333333 - 0.0454545 + 0.006060 
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D2 = 0.0740342 + 0.0149190 - 0.00118583 

- 0.0045481 . 0.00008286 + 0.0012013 L^L2 

D3 = 0.08461538 + 0.02435846 - 0.00293040 

D4 = 0.4174603 

Wilks made the assumption that a basically linear relationship 
existed between the boundary-layer thickness and the momentum 
thickness. This assumes that D2 does not vary significantly 
with X. A basically linear relationship between the mixed- 
momentum thickness and the boundary-layer thickness was also 
assumed . 

With the introduction of a dimensionless parameter 
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the foregoing equations can now be written 
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Combining these three equations yields 
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In addition to the above relationships, the previously derived 
equation for the shape factor L 2 is also important, viz.. 




The following characteristics must be considered for the 

perfectly conical nozzle. 

r(X) = r (1 - X) 
o 

By continuity 
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Conservation of angular momentum requires that 

Vo*rQ "" V(X) ‘r^* (1 - X) = A 

If V = K*U^ 
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APPENDIX B 



DERIVATION OF PRESENT THEORY 



The concept of a free vortex in the inviscid core requires 
that the tangential velocity be inversely proportional to the 
radius. In the case where boundary layer growth is large, the 
boundary layer thickness can no longer be considered thin in 
comparison to the radius. VJith these considerations, the 
boundary conditions for the momentum equation become 
V V(x) = A /(r-6) as z -> <» 

This new definition of V (x) has no effect upon the development 
of the u momentum-integral equation. However, there is an 
effect in the v roomentxim-integral equation. The development 
is the same as described in Appendix A up to the point 
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Thus the momentum-integral equation becomes 




dr u* V dr u« V 

dx r dx r-6 





69 



Simplif ing 
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Introducing the velocity profiles as in Appendix A, the momentum- 
integral equations and the appropriate boundary conditions 
become 



hk f(f-l)-an) + (f- 1 ) 

+ _L_-v 2.^.6- (l-g2)dn = - v (0) 

c • r uA Jo ° 



■dn 



with the boundary conditions 



(i) 



f = 0 



/• -X , U dU 

(XI) - y? f’ = c'd“ 



dr 

c*r d 



on n = 0 



and 



(iii) 

(iv) 

(V) 

1 d 



f = 1 

f ' = 0 

f" = 0 



> on n = 1 



Ido U*V 1 /d6 6 dr\ 

cr^ dx Jo r-6 c \dX r dx; j 



cr 

= - ^9'(0) 

with the boundary conditions 
(i) g = 0 

(ii) . g" = 0 



on n = 0 
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I 



« 









4 



(iii) g = 1 

(iv) g' = 0 

(v) g" = 0 



►on n = 0 



Again, it was assumed that a basically linear relationship 
existed between the boundary-layer thickness and both the 
momentum and mixed -momentum thicknesses. Defining the 
coefficients Al, A2, A3, and A4 as in Appendix A, and intro- 
ducing a new term 



A5 = 



fl 



f dn 



f (l-g)dn 



1 - D1 
D3 



the momentxam-integral equations then became 



r ' U ' 6 ' 

~ + U“ (2 + Al) + ^ 



A2 



r* /VV 
r *lu) 



•A3 



2r' U' V . 6’ . 1 f6' 

r U V "^6 r/6 - 1 \6 




A4 

G 



The second equation can be rewritten 




/(I + 



A5 



) 



Thus, the shape factor can be obtained from 

A2 = gF^ . (2 + Al) + —• (1 - ^ A3) + — 1 
Lu r 6 J 



and the nev; values for G and L 2 can be calculated from the 
value of G', where G' is now 
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I 





and 




As the boundary layer grows, the axial velocity as defined 
in APPENDIX A for a cone will no longer be valid. The increas- 
ing boundary-layer thickness will cause the flow in the axial 
direction to accelerate more than would be caused by the 
convergence of the cone alone. The concept of displacement 
thickness would take into consideration this change in the 
effective radius of the cone. Continuity now yields 



“o'['^o ■ = U(X)-[^r^-a-X) - ^-TT 



thus 



and 



U 



U 



U(X) = 



1 - 






] 



(1-X) - 









(1-X) - 



6l^(X) 



With the new definition of V (X) , the conservation of angular 
momentum yields 
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V • r - 6(0) = V(X)- r ‘(l-x) - 6 (X) 

o o ° 
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